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When NMR lines overlap and at least one of them is affected 
by radiation damping, the resonance line shapes of all lines are 
no longer Lorentzian. We report the appearance of narrow 
signal distortions, which resemble hole-burnt spectra. This new 
experimental phenomenon facilitates the detection of tiny sig- 
nals hidden below the main resonance. Theoretical analysis 
based on modified Maxwell-Bloch equations shows that the 



presence of strong transverse magnetization creates a feed- 
back through the coil, which influences the magnetization of 
all spins with overlapping resonance lines. In the time domain 
this leads to cross-precession terms between magnetization 
densities, which ultimately cause non-linear behavior. Numeri- 
cal simulations corroborate this interpretation. 



1. Introduction 

Radiation damping is an effect known since the early days of 
NMR. [1] This non-linear phenomenon results from the interac- 
tion of the magnetization with itself, mediated through the 
tightly coupled detection coil. Radiation damping is nowadays 
routinely encountered in high-resolution liquid-state NMR 
spectroscopy due to the increase of magnetization achieved 
by high static magnetic fields and/or hyperpolarized species 
and by the improvements of the electronic detection circuit, in 
particular by cold probes, which exhibit huge quality factors. 
Radiation damping is often seen as a spurious and undesired 
effect which broadens resonance lines for small flip-angle exci- 
tation pulses, [1] distorts the line shape for larger flip angle 
pulses, [2] and even prevents the safe and reliable inversion of 
longitudinal magnetization. [3] As a consequence several NMR 
methods, such as WATERGATE/ 41 water flip back, [5] presatura- 
tion, [6] jump-and-return [7] or electronic devices for controlling 
radiation damping [8] have been developed either for taking 
benefit of the fast transverse decay rates of magnetization ex- 
periencing radiation damping or for suppressing its effects for 
intense signals from the acquired spectra. More recently the 
situation has changed. Instead of being considered as an un- 
pleasant effect, that has to be avoided, its non-linear feature is 
used for experimental studies, either for exploring new physi- 
cal phenomena such as chaotic behavior, [9] potentially en- 
hanced by electronic circuits/ 101 or multiple spontaneous maser 
emissions [11] or for evaluating alternative detection schemes 
such as radiation-damping amplified signal detection [12] or nu- 
clear-spin noise. [13_15] The present study belongs to this last 
class. Indeed we have experimentally observed that the line 
shapes of spin systems experiencing radiation damping and 
that have close resonance frequencies, are strongly affected 
and strongly differ from the usual bell-shaped curves observed 
for small excitation pulses. Indeed, they seem to be affected 
by a hole-burning phenomenon. We show that this effect re- 
sults from additional cross-precession terms present in the 
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spin dynamics, which couple the evolutions of different mag- 
netization densities through the induction they create inside 
the coil. Thanks to these non-linear interactions, we show that 
tiny signals resonating within the line width of large ones 
become easily detectable, since they introduce strong narrow 
distortions into the line shape, which are amplified or look like 
"holes" in the main resonance line instead of the expected 
small positive peaks or shoulders. 
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2. Theory 

2.1. Radiation Damping Field 

We consider the simple resonant electronic detection circuit 
depicted in Figure 1 and follow the derivation of Vlassenbroek 
et al. for computing the radiation damping field. [16] We assume 
that the detection coil is perpendicular to the static magnetic 
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Figure 1. Classical probe with a coil of inductance L, resistance R, and with 
tuning capacitor C. 

field B 0 , that we assume to be perfectly homogeneous. As a 
result, the transverse component of the precessing nuclear 
magnetization Ji x (t) creates an electromagnetic force (emf) 
V s {t) in this detection coil, which can be calculated thanks to 
the reciprocity principle given by Equation (1): 



sequently consider only the steady-state AC regime for the so- 
lution of Equation (4): 

Uw + r + -^jl(t) =Z 0 (co)l(t) = -y/wLr^Jf^t) (4) 

with Z 0 (oj) = |Z 0 (a;)|e #(ft>) , the complex impedance [Eq (5)]: 
Lco LC 



l^oWl = + a2 lc = arctanA LC 



A LC is a function of the offset between the Larmor frequency 
and the electronic circuit resonance frequency Eq. (6): 



CO LC CO 



CO 



(6) 



Finally the magnetic field B q (f) created by the precessing 
transverse magnetization is given by Equation (7): 



(7) 



Using u x , u y , u z as the unit vectors of the respective coordi- 
nate axes, the equations describing the evolution of the nucle- 
ar magnetization density Jl are obtained by standard theo- 
retical approaches. Considering a single nuclear spin 1/2 
system, which can be described by the Bloch equations, we 
have Equation (8): 



V s (t) = - y /fi 0 rjLr- rt Jt x (t) 



m 



where ft 0 is the the magnetic permeability of free space, rj the 
filling factor, L the inductance of the coil, and "V the sample 
volume. This emf creates a current /(f) in the idealized detec- 
tion circuit [Eq. (2)]: 



d 2 



+ /?^/(t) + ^/(r) 



dt 



By introducing the characteristic parameters co lc = '\/VlC l 
the resonant frequency of the detection circuit, and Q = Lco lc /R, 
the quality factor of the coil, the equation describing the flow- 
ing current as a function of the precessing transverse magneti- 
zation is obtained [Eq. (3)]: 



d? /(f) + ^S /(f) + ^ /(f) 



L dt 2 ' 



This differential equation corresponds to that of a damped 
harmonic oscillator driven by an external source. Transient 
phenomena are present for periods of the order of a few 
Q/co lo which typically represents a few tens of microseconds 
for liquid-state NMR probes with the highest Q. This is most of 
the time a short duration compared to the typical evolution 
time of magnetization (hundreds of milliseconds). We can con- 



dt 



(f) =yJt{t) x BotC+B^f) 



z(0 



(8) 



where y is the magnetogyric ratio, and T } and T 2 are the longi- 
tudinal and transverse self-relaxation times. Equations (7) and 
(8) are the Maxwell-Bloch equations. Since B q (t) is dependent 
on the transverse magnetization, the evolution of the magneti- 
zation, as described by Equation (8) is no longer linear. These 
equations can be numerically integrated to explore the effects 
of distributions of Larmor resonance frequencies or the effects 
of a mis-tuned electronic circuit. [3,9,16] Several studies have 
been reported in these cases and in particular an analytical so- 
lution for negligible longitudinal relaxation and perfect 
tuning (A LC = 0) even exists whatever the initial flip-angle exci- 
tation pulse. [2] 



2.2. Radiation Damping in the Presence of Two Resonance 
Lines 

Herein we consider the existence of two different species of 
magnetization {Jt° and M b ) of different Larmor resonance fre- 
quencies {yC\-cf)B 0 and o b )B 0 , where o° and d° are the 
shielding constants of Ji° and Jt h , respectively). The principles 
used in the preceding derivation are still valid, and we can 
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deduce that the magnetic field B-,(f) is now given by Equa- 
tion (9): 1 



B,(f) 



fi 0 rjL ( d 



(9) 



Ji b Larmor frequency, o^Bq. As usual for a cross-preces- 
sion term between transverse components, it will be complete- 
ly negligible, that is, non-secular, if the difference of resonance 
frequencies \ y{o a —o b )B 0 ) \ is much larger than the typical radia- 
tion damping characteristic rate Ap D induced by Ji b [Eq. (13)]: 



and the evolution of the magnetization, still considering inde- 
pendent spins 1/2, is obtained by inserting the new expression 
of B^f) into Equation (8), to yield Equations (10) and (11): 



2\Z n (a) ir )\ 



(13) 



dt 



(f) - <f) Ji (f) x 



fi 0 rjL ( d 



T? 



(10) 



^ Jf (f) =y(l -a b ) .A (f) x 



'ZJw) \dt 



(11) 



7? 



where straightforward extensions of the notations were used. 
These six differential equations [Eqs. (10) and (11)] are all cou- 
pled together through the radiation damping field [B-,(f), 
Eq. (9)], although we have considered independent spins. The 
evolution in the rotating frame is obtained by considering only 
one of the two rotating components of B-,(f). Finally in a frame 
rotating at co 0 = yB 0 , this coupling between the evolutions of 
the transverse magnetization densities of Ji° and Ji b is readily 
illustrated by the evolution of one rotating component 
[Eq.(12)]: 



This corresponds to the maximal transverse decay 
rate [contribution from Jt b to the B q (f) field], obtained 
for perfectly tuning A LC = 0, and very small flip angles 
[Eq. (12)]. [12] Ap D can be estimated from the full width 
at half maximum (FWHM in Hz) of the Ji b resonance 
after a small flip excitation pulse, since (X b RD + 1 /T 2 b ) = 
\ 1 JtFWHM. In summary, this cross-precession term has to 
^x(0 j be considered if the two resonance lines of Jl° and 
M b are overlapping or are close to this condition. 

Although this semi-quantitative description provides 
the experimental conditions where cross-precession 
between two independent resonances due to radia- 
tion damping can be observed, the analytical prediction of the 
resonance line-shapes is prevented by the non-linear charac- 
ters of Equations (10) and (11). Numerical simulations, however, 
can provide a detailed insight into the resonance line shapes. 



2.3. Extension of the Theoretical Framework 

Several comments lead us to extend the previous derivation to 
take into account more complex spin systems or other contri- 
butions neglected in a first approximation: 



fi 0 t]LY 2 B 0 
2Z 0 (w) • 



t - 



<(0 
7? 



(12) 



In this equation we recognize the Zeeman term, the classical 
radiation damping term, the cross-precession term due to 
feedback field and finally the transverse relaxation contribu- 
tion. In summary, through the coupling to the detection coil 
cross-precession terms between the two magnetization densi- 
ties appear. 

Z 0 (co) being a slowly varying function of co around the 
Larmor resonance frequencies, since it depends on the half- 
width of the tuned electronic circuit at 3dB, [cd ic /(2Q)], can be 
considered as constant and the two contributions to the B q (f) 
field, that is, the cross-precession terms, can be compared di- 
rectly. Considering the M° magnetization, the cross-precession 
term due to radiation damping resulting from the transverse 
magnetization of Jt b can be seen as an off-resonance rf excita- 
tion whose amplitude is proportional to Jt b x and whose fre- 
quency is equal, to a first approximation (see below), to the 

7 Note that contrary to ref. [3] we explicitly consider the effect of different reso- 
nance frequencies in this derivation, which is needed due to the reduced am- 
plitude of B, (f). 



The extension from two spins to n species is obtained by 
adding their contributions to the B-, (t) field [Eq. (9)] and by 
modifying the expression of Bloch equations [Eq. (10)] ac- 
cordingly. 

As a first approach, for a simple 1D acquisition spectrum, a 
scalar-coupled m-multiplet can be described as a superposi- 
tion of m fictitious spins 1/2 with appropriate intensities, 
the evolution of which follows the previous description. 
According to the tuning frequency of the coil a> lo the angle 
ipico) between the radiation damping field and the magneti- 
zation varies as well as the coupling between the magneti- 
zation and the coil [Eq. (5)]. [14,16] This gives rise to a variation 
of the apparent radiation damping characteristic rate and 
to a "cavity pulling" effect or more rigorously to a frequency 
pushing. [14] The latter leads to a variation of the central fre- 
quency of the line affected by radiation damping. [17,18] This 
frequency shift results from the imaginary component of 
Z 0 ((o), and is thus on the order of sin^MA RD , that is, up to a 
few hundreds rads" 1 for experiments performed with a 
cold probe. 

A strong magnetization creates an average dipolar field, 
which superimposes on the static magnetic field and is 
thus experienced by all spins. [19] The amplitude of this aver- 
age dipolar field is dependent on the z component of the 
magnetization, on the sample shape, and on a factor of 3/2 
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or 1 depending on the nature of the magnetization which 
creates the field ("like" or "unlike" spins [20] ). These contribu- 
tions should be added to the Bloch equations, and are typi- 
cally time-dependent during the acquisition (a variation of 
the longitudinal magnetization potentially counterbalanced 
by the B 0 feedback field controlled by the lock signal). The 
order of magnitude of these contributions {~ft 0 yJ? 0 ) is less 
than 10 rads" 1 for a bulk solvent in a high magnetic field at 
room temperature. It can consequently easily be masked by 
cavity pulling (~A lc rjQ/2xju 0 yJ? 0 ). Herein, we neglect these 
average dipolar field contributions in the numerical simula- 
tions, in part because we used small flip-angle excitation 
pulses, so that Jt^M^ is nearly constant. 
• Finally transverse and longitudinal dipolar cross-relaxations 
(and dynamic frequency shifts) between the different mag- 
netization species can be included. These contributions are 
typically one tenth of the intermolecular dipolar self-relaxa- 
tion contributions, that is, on the order of 0.1 Hz and conse- 
quently they have not been included in the numerical simu- 
lations, here. For systems of equivalent spins, in particular 
methyl groups, with long correlation times, cross-correlation 
effects [21] may be larger and also need to be considered for 
further refinement. 



3. Results and Discussion 

Figure 2 illustrates the types of unexpected line-shapes ob- 
served when two different spin species with similar chemical 
shifts are present and at least one of them is significantly influ- 
enced by radiation damping. In Figure 2a, the two species 
have identical intensities since they correspond to the methyl 
doublet of isopropyl acetate. For two overlapping, broadened 
Lorentzian lines one would usually expect symmetrical broad- 
ening of each line, eventually leading to their coalescence and 
to the disappearance of the doublet structure. Instead we ob- 
serve that each line is strongly asymmetric and a steep de- 
crease in intensity occurs between them. This appearance of a 
dip or hole in the line shape is even more obvious in Fig- 
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Figure 2. Example of non-linear response induced by radiation damping 
when two lines overlap. These spectra were acquired in conditions where 
cavity pulling vanishes bp{a>) = 0] with small flip-angle excitation pulses, 
a) Case of the doublet of isopropyl acetate. This spectrum was acquired at 
600 MHz with a cold probe. The solution contained 50% vol/vol of isopropyl 
acetate, b) Case of a mixture of acetone and acetonitrile (40:5) with a con- 
centration of 80% in volume of acetone (lock solvent: [D 6 ]DMSO). This spec- 
trum was acquired at 700 MHz with a cold probe. Experimental signals 
marked by an asterisk are 13 C satellites. 



ure2b, where we used a mixture of acetone and acetonitrile 
(40:5). Due to radiation damping broadening, the acetone res- 
onance overlaps the acetonitrile one. Instead of observing a 
small peak or a shoulder at the acetonitrile chemical shift, 
which would be difficult to identify, one can notice an easily 
detectable dip, the magnitude of which is not directly correlat- 
ed to the ratio of concentrations. 

In order to check the consistency of the proposed interpre- 
tation, a series of 1D spectra of a mixture of acetone and ace- 
tonitrile (40:5, total acetone concentration 80% in volume) 
was acquired with small flip angle excitation pulses at different 
setting of the resonance frequency of the rf-circuit (co lc ). The 
results of these experiments are shown Figure 3 a. Numerical 
simulations corresponding to this series of experiments are re- 
ported in Figure 3 b. They were obtained by numerically inte- 
grating the system of differential equations [Eqs. (10) and (11) 
expressed in the rotating frame] using parameters [Ji Ql Q, rj) 
deduced from the experimental conditions. 

Several items can be seen from Figure (3): 

• The comparison between the experimental and simulated 
data reveals near-perfect agreement. This clearly substanti- 
ates the theoretical derivation and proves that the non- 
linear coupling between the two magnetization densities 
mediated by the detection coil is the reason for observing 
such unexpected line shapes. 

• The peak shape and position of the acetone signal are 
strongly dependent on the tuning conditions, varying over 
30 Hz. This effect results from cavity pulling. [17,18] As a func- 
tion of the amount of this contribution the resonance lines 
of acetone and acetonitrile are more or less overlapping. By 
contrast, the smaller concentration of acetonitrile reduces 
its magnetization and thus the amplitude of the radiation 
damping field [Eq. (13)]. As a consequence the position of 
this line is practically unaffected by cavity pulling. 

• In addition to the variation of the acetone peak's resonance 
frequency, its line-width is dependent on the tuning offset. 
It is broadest when the cavity pulling effect vanishes (red 
curves in Figure 3). At this condition we observed a probe- 
dependent offset from the point where the sensitivity of 
the detection circuit was previously found to be highest 
(the spin noise tuning optimum, SNTO). [22,23] This offset re- 
sults from the fact that, while the feed-back field is gov- 
erned by the flowing current inside the coil, the SNTO cor- 
responds to a voltage optimum. The high impedance of the 
preamplifier then causes the difference between the two 
resonance situations. While these conditions allow highly 
sensitive detection of signals not or only moderately affect- 
ed by radiation damping, an intense signal of long intrinsic 
transverse self-relaxation time T 2 is reduced in sensitivity 
due to the extremely strong radiation damping. These two 
observations are in agreement with the fact that the in- 
duced change of Aco LC affects the radiation damping rate: 
in a first approximation (small offsets), it broadens the line 
as A RD cos^M and shifts it by cavity pulling as X RD s\nip(co), 
with A RD taken at co LC = yB 0 . 
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The system of non-linear equations [Eqs. (10) and (11)] is 
valid whatever the flip angle of the excitation pulse. In the 
present context, to avoid unreliable behavior of the electronic 
circuits (non-linearity of the excitation line due to mismatch, 
saturation of the preamplifier or of the ADC) we have limited 
the experiments to small flip-angle excitation pulses. This 
domain is also physically relevant because it corresponds to 
the one where linear-response theory is usually valid. [24] 
Beyond this domain, the resonance line shapes can be calculat- 
ed using Equations (10) and (11) for any flip-angle excitation 
pulses. Beside the severe deviations from the Lorentzian 
shapes due to radiation damping as already reported for a 
single line, [2] this investigation reveals the appearance of the 
same peculiarities (amplification, asymmetric shape, dip in the 
line shape) due to cross-precession, which superimpose on re- 
sonance lines. 

The agreement between the experimental and simulated 
curves allows the evaluation of this radiation-damping cross- 
precession term for spectroscopic applications. Indeed, the 
strong distortions observed on the line shape, when resonan- 
ces overlapped, in which at least one of them is affected by ra- 
diation damping, provide routes for improved detection. In 
Figures 4 a,a' we illustrate this by the results obtained by nu- 
merical simulations of two species with a ratio of concentra- 
tions of 100:1. On the one hand, when one considers the 
system with a linear response, that is with a probe of small 
quality factor Q (Figure 4a'), one can notice that the smallest 
peak is not detected. On the other hand, for a probe with a 
large Q, that is when radiation damping is present the non- 
linear response of the system allows simple detection of the 



Figure 3. Variation of the resonance line shape of a mixture of acetone and 
acetonitrile (40:5) with a concentration of 80% in volume of acetone as a 
function of the resonance frequency of the electronic circuit a> LC . a) Experi- 
mental spectra obtained on a 500 MHz Bruker Avance III spectrometer 
equipped with a TXI cryoprobe. Different colors represent different offsets 
A LC relative to the condition where oj lc matches the Larmor frequency 
(shown in red). For taking into account the effect of different detection sen- 
sitivities due to offset tuning, all spectra were normalized by using the 
DMSO signal whose concentration was sufficiently low so as to not being af- 
fected by radiation damping. Signals marked by asterisks are 13 C satellites, 
b) Numerical simulations obtained by integrating Equations (10) and (11) 
using parameters {Jt 0l Q, rj, A LC ) corresponding to the experimental cases. 



The shape of the acetonitrile resonance line appears to be 
strongly dependent on that of acetone. Indeed when the 
latter is relatively narrow so they do not overlap, the aceto- 
nitrile resonance line is of Lorentzian shape and may be 
amplified, as easily observed on the 13 C satellites. Increasing 
the overlap and thus the contribution from radiation damp- 
ing to the broadening of the acetone line, the acetonitrile 
resonance changes from a peak superimposed on the ace- 
tone resonance to a dip or hole within the broad peak. 
Based on the numerical simulations under "hole condi- 
tions," the exact resonance frequency of the minor compo- 
nent appears not to be equal to that of the minimum but it 
is slightly shifted in the direction opposite to the major 
peak's maximum. 
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Figure 4. Numerical simulations of different conditions allowing one to see 
the influence of the existence of radiation-damping cross-precession term, 
a) and a') Case of two signals with a proton ratio of 100:1. Without radiation 
damping (a') the signal corresponding to the smaller signal can hardly be 
detected, while with radiation damping (a) this becomes possible due to the 
appearance of a dip in the line, b) and b') In case of radiation damping 
acting on a symmetrical doublet, a hole appears between the two lines (b). 
In the simulations without radiation damping (b') the two lines are broad- 
ened to give approximately the same line width as in (b) by using small 
values for T 2 . This leads to the overlap of two Lorentzian lines and the disap- 
pearance of the doublet structure. 
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smaller peak since a narrow distortion in the main resonance is 
clearly visible. 

Another example of the unique features of this radiation- 
damping cross-precession phenomenon is provided in Figur- 
es 4 0,0', where two resonances of identical intensity are con- 
sidered. There, the transverse relaxation rates are adapted such 
that they exhibit the same line-widths for non-overlapping 
lines. In Figure 4b' a low-Q probe is assumed so that no radia- 
tion damping broadening or cross-precession is present. As 
one can notice, a large unstructured peak results. By contrast, 
in Figure 4 b, a high-Q probe is assumed, that is, cross-preces- 
sion is now included. The line shape is strongly affected since 
a hole in the middle is observed. This provides a way of clearly 
assigning the existence of two lines, which would be rather 
difficult in the first case. 

4. Conclusions 

Herein, we showed that the resonance line shapes of spin sys- 
tems experiencing radiation damping can exhibit severe distor- 
tions when resonances overlap. In particular, instead of the ex- 
pected superposition of positive peaks when two lines overlap, 
one can observe a significant dip within the line shape. This 
previously unknown effect is explained by the fact, that 
through the detection coil, the precessing magnetization den- 
sities create a feedback field, which is experienced by all reso- 
nances. As a result cross-precession terms dependent on the 
different magnetization amplitudes have to be included in the 
Bloch description of the spin dynamics. These terms become 
non-secular when the resonance frequencies are sufficiently 
different. The deduced non-linear differential equations are 
able to represent all the experimentally observed behaviors. 
The discovery of this unexpected effect and its theoretical in- 
terpretation obviously raises the question of its practical appli- 
cation. First, the incidence of radiation-damping-broadened 
lines is expected to increase in the near future due to the avail- 
ability of cold probes, which have very large quality factors Q 
and due to the increase of the magnetization achievable by 
high magnetic fields and/or by hyperpolarized species. Second, 
as illustrated by the simulations shown in Figure 4, the non- 
linear response induced by this cross-precession effect allows 
the improvement of detection capabilities through spectral 
features (existence of a doublet, tiny signals, etc.) which would 
have been nearly undetectable in the linear regime. 

Experimental Section 

All experiments reported herein were performed at room tempera- 
ture using Bruker Avance spectrometers running at 500, 600 or 
700 MHz equipped with cryoprobes. Similar spectra were also ob- 
tained with classical room-temperature probes. Samples were com- 
posed of mixtures of different solvents of analytical grade. Deuter- 
ated DMSO was also added for field frequency locking and to con- 
trol total spin densities. ID NMR spectra were acquired using small 
flip-angle excitation pulses (< 15°). 

The spectral simulations were obtained, using Scilab software, by 
Fourier transforming the simulated time-domain signal amplitudes, 



which were computed by numerically integrating the system of 
Equations (10) and (11) expressed in the rotating frame. In order to 
have simulated data as closely comparable as possible to the ex- 
perimental ones, the different parameters were adjusted to the re- 
spective experiments. In particular the proton concentrations 
needed for computing the total magnetization were those of the 
solution, the quality factor Q was set equal to 550 for the 500 MHz 
spectrometer in agreement with cavity-pulling experiments/ 181 and 
direct Q measurements performed with a spectrum analyzer. From 
the experimental line-width a value of X RD of 62 Hz on the 500 MHz 
and from the Q value and the magnetization, a filling factor rj of 
0.03 was deduced. The latter value is in agreement with previous 
reports of filling factors for saddle coil probes. [25,26] 
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